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PREFACE 


This  report  is  a  presentation  and  analysis  of  the  electrical  resistivity  and  magne¬ 
tometer  airborne  surveys  flown  in  September  and  October  of  1975  in  coordination 
with  the  proposed  Dickey-Lincoln  School  Dam  Project  in  northern  Maine.  It  was 
prepared  by  Paul  V.  Sellmann,  Geologist,  of  the  Northern  Engineering  Research 
Branch,  Experimental  Engineering  Division,  and  by  Steven  A.  Arcone,  Geo¬ 
physicist,  and  Allan  J.  Delaney,  Physical  Sciences  Technician,  of  the  Physical 
Sciences  Branch,  Research  Division,  U.S.  Army  Cold  Regions  Research  and  Engi¬ 
neering  Laboratory.  The  Foundations  and  Materials  Branch  of  the  U.S.  Army 
Engineer  Division,  New  England  funded  the  research.  Previous  CRREL  participation 
in  this  project  includes  investigation  of  the  surficial  geology  of  the  area  and  tectonic 
activity  of  the  region  (McKim  and  Merry  1975). 

The  airborne  surveys  were  flown  by  Barringer  Research  Ltd.  of  Toronto, 

Canada,  under  contract  to  the  New  England  Division  of  the  U.S.  Army  Corps  of 
Engineers.  The  multi-elevation  surveys  and  ground  studies  were  supported  by 
DA  Project  4A76271 9AT24,  Design,  Construction  and  Operation  Technology  for 
Cold  Regions,  Task  A2,  Soils  and  Foundation  Technology  for  Cold  Regions,  Work 
Unit  003,  Electromagnetic  Methods  for  Subsurface  Exploration  (OCR  1 .07-CARDS 
114).  Upon  completion  of  the  surveys  all  data,  maps,  flight  records,  photography, 
and  daily  logs  were  sent  to  CRREL  where  further  processing,  evaluation,  and  geo¬ 
logic  interpretations  and  correlations  were  performed. 

Pieter  Hoekstra  and  Frank  jagodits  technically  reviewed  the  manuscript. 

The  contents  of  this  report  are  not  to  be  used  for  advertising  or  promotional 
purposes.  Citation  of  trade  names  does  not  const!  iute  an  official  endorsement  or 
approval  of  the  use  of  such  commercial  products. 
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SUMMARY 


Simultaneous  airborne  magnetometer  and  VLF  resistivity  surveys  were  con¬ 
ducted  in  northern  Maine  to  locate  rock  types  suitable  for  construction  purposes. 
This  was  done  for  the  New  England  Division  of  the  U.S.  Army  Corps  of  Engineers 
in  response  to  anticipated  requirements  for  the  proposed  Dickey-Lincoln  School 
Dam  Project.  The  surveys  covered  approximately  920  km2 .  In  addition,  a  small, 
special  test  area  of  known  bedrock  geology  was  surveyed  to  establish  specific 
ground  correlation  and  to  examine  the  effects  of  flight  altitude  upon  resolution. 

Over  the  special  test  area  both  the  contoured  magnetometer  and  resistivity 
dat?  differentiated  the  known  igneous  geology  from  the  surrounding  slate.  Re¬ 
sistivity  data  obtained  by  ground  VLF  methods  correlated  well  with  the  airborne 
resistivity  results.  The  magnetic  data  indicated  a  normally  expected  degree  of 
magnetic  mineralization  associated  with  the  granodiorite  and  syenite  in  the  test 
area,  which  can  contain  magnetite  as  an  accessory  mineral  (Boone  1962).  The 
multi-elevation  surveys  were  flown  at  mean  flight  altitudes  of  1 50  m  and  300  m. 
The  general  distribution  of  anomalies  at  the  higher  altitudes  remained  unchanged, 
although  detail  was  decreased.  These  results  implied  that  local  changes  in  altitude 
occurring  during  this  survey  would  not  cause  any  significant  loss  in  resolution.  In 
general,  the  test  area  results  ensured  the  reliability  of  the  overall  survey  and  es¬ 
tablished  a  basis  for  airborne  identification  of  igneous  rocks  for  this  study. 

The  data  from  the  entire  survey  indicated  the  following  results: 

1 )  No  other  resistivity  anomalies  occurred  with  values  as  high  as  those  associated 
with  the  igneous  rocks  of  the  test  area. 

2)  No  magnetic  anomalies  were  found  outside  the  test  sector. 

3)  The  more  resistive  areas  outside  the  test  sector  corresponded  with  a  cyclically 
bedded  gray  slate  and  sandstone  unit  in  the  central  part  of  the  survey  area. 

4)  The  resistivity  data  corresponded  well  with  mapped  reconnaissance  geology 
of  the  region  (Boudette  et  al.  1966). 

From  these  results  it  was  concluded  that  no  new  sources  of  igneous  rocks  are 
apparent.  Rocks  most  suitable  for  construction  purposes  near  the  dam  site  may 
occur  in  the  central  part  of  the  survey  area,  which  is  associated  with  areas  of  high 
resistivity. 
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AIRBORNE  RESISTIVITY  AND  MAGNETOMETER 
SURVEY  IN  NORTHERN  MAINE  FOR  OBTAINING 
INFORMATION  ON  BEDROCK  GEOLOGY 


P.V.  Sellmann,  S.A.  Arcone  and  A.J.  Delaney 

INTRODUCTION 

Knowledge  of  the  distribution  and  quality  of  bed¬ 
rock  within  the  area  of  the  proposed  Dickey-Lincoln 
School  Dam  Project  is  important  for  development  of 
construction  plans  and  environmental  impact  statements. 
At  the  conception  of  this  study  known  sources  of 
rock  suitable  for  construction  purposes  were  of  limited 
extent,  since  the  most  common  rock  type  in  the  area 
is  slate.  The  most  suitable  rocks  are  part  of  a  remotely 
situated  intrusive  body  more  than  19  kilometers 
from  the  proposed  dam  site.  Background  geologic 
information  is  available  from  a  regional  reconnaissance 
study  and  a  detailed  investigation  of  the  known  intru¬ 
sive  complex  (Boudette  etal.  1966,  Boone  1962). 

Based  on  these  studies  the  New  England  Division  of 
the  Corps  of  Engineers  (NED)  felt  that  additional  data 
were  required  to  determine  if  unknown  sources  of 
rock  suitable  for  construction  purposes  could  be  locat¬ 
ed  nearer  the  project  sites.  As  a  result,  CRREL  was  ap¬ 
proached  to  determine  if  resistivity  techniques  could 
be  used  as  an  aid  in  obtaining  this  information. 

In  March  1975  a  preliminary  study  was  made  to  deter¬ 
mine  if  resistivity  contrasts  of  known  rock  types  in  the 
region  were  great  enough  to  justify  an  extensive  ground  or 
airborne  survey  using  the  radiowave  resistivity  tech¬ 
nique  at  VLF  (very  low  frequency)  which  is  the  fre¬ 
quency  range  best  suited  for  a  bedrock  study.  It  in¬ 
dicated  that  contrasts  in  the  electrical  properties  be¬ 
tween  the  intrusive  rocks  and  other  rock  types  in  the 
area  were  great  enough  to  permit  this  method  to  be 
considered  for  the  proposed  purpose,  as  well  as  for 
providing  additional  general  information  on  bedrock 
properties  (Sellmann  et  a!.  1974, 1975).  As  a  result 


of  this  work  CRREL  recommended  that  additional 
field  investigations  were  warranted,  and  both  an  air¬ 
borne  resistivity  and  magnetometer  survey  were  con¬ 
tracted.  The  surveys  were  used  to  complement 
each  other,  providing  both  electrical  and  magnetic 
property  data  for  the  rock  types  as  well  as  a  means 
of  detecting  magnetic  mineralization. 

The  extent  of  airborne  survey  coverage  was  jointly 
determined  by  NED  and  CRREL  project  personnel 
(Fig.  1 ).  Because  of  the  large  size  and  irregular  shape 
of  the  total  reservoir  impoundment  limits,  complete 
airborne  survey  coverage  was  not  practical  or  con¬ 
sidered  necessary.  With  the  emphasis  of  this  project 
placed  on  obtaining  bedrock  data,  the  more  central 
part  of  the  area  near  the  proposed  major  construction 
activity  was  covered  in  detail.  For  control  purposes 
these  limits  were  also  adjusted  to  include  all  pre¬ 
viously  mapped  major  rock  types.  Other  factors 
influenced  shape,  e.g.  the  minimum  practical  length 
of  flight  lines,  which  is  approximately  16  km,  and 
the  spacing  of  flight  lines.  A  flight  line  spacing  of 
0.4  km  was  used  for  this  study,  in  contrast  to  a  tenth 
of  a  mile  used  on  some  previous  surveys  (Hoekstra 
etal.  1974). 

In  areas  outside  the  main  survey,  pairs  of  magne¬ 
tometer  profiles  were  flown  along  the  drainage  net¬ 
works.  Only  magnetometer  data  were  obtained,  since 
these  can  be  acquired  independently  of  flight  line  ori¬ 
entation,  while  resistivity  lines  must  be  flown  at  a 
fixed  orientation  for  maximum  coupling  with  the  re¬ 
mote  VLF  transmitter. 

During  the  first  week  of  the  airborne  survey  ex¬ 
tensive  ground  measurements  were  performed  in  a 
special  test  sector  of  known  bedrock  geology,  shown 
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Figure  1.  Index  map  showing  the  extent  of  airborne  coverage.  The  survey  area  was  flown  with  both  the  re¬ 
sistivity  and  magnetometer  systems.  Pairs  of  magnetometer  lines  were  flown  along  the  major  drainage  net¬ 
works. 
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a.  Ground  unit  (Geonlcs  EM16R)  for  measuring  resis¬ 
tivity  at  VLF by  using  the  surface  Impedance  method. 


b.  STOL  aircraft  with  E-PHASE  antennas  In  the  nose  cone  and  a  magnetometer  sensor  extended  to  the  rear 
(photograph  courtesy  of  Barringer  Research  Ltd.). 

Figure  2.  Systems  for  measuring  resistivity  at  VLF. 


in  Figure  1 .  In  addition  two  separate  surveys  in  this 
sector  were  conducted  at  mean  altitudes  of  1  SO  and 
300  m.  These  special  studies  were  made  to  verify  the 
quality  of  the  entire  airborne  survey,  with  respect  to 
both  rock  type  differentiation  and  the  effect  of  flight 
altitude  upon  resolution  of  resistivity  anomalies. 


MEASUREMENT  TECHNIQUES  EMPLOYED 
Ground 

A  commercially  available  Geonics  EMI  6R  ground 
unit  for  measuring  resistivity,  illustrated  in  Figure  2a 
and  on  the  cover  of  this  report,  wai  used  for  monitor¬ 
ing  the  VLF  transmitter  (NAA,  17.4'  kHz)  located  near 
Cutler,  Maine.  The  instrument  is  calibrated  directly  in 
apparent  resistivity  and  phase  for  the  particular  fre¬ 
quency  desired.  The  measurement  is  performed  by 
determining  the  ratio  of  Ex  (the  horizontal  electric 
field),  measured  between  two  probes  spaced  at  10  m 
and  oriented  in  a  radial  direction  towards  the  trans¬ 


mitter,  and  Wy  (the  horizontal  magnetic  field),  meas¬ 
ured  with  a  small  coil  near  the  ground  surface 
orthogonally  to  fx.  Since  no  current  need  flow  from 
the  ground  to  the  probes  the  system  is  virtually  free 
of  contact  resistance  problems. 

The  coil  is  used  to  find  the  direction  to  the  station, 
and  the  correct  apparent  resistivity  values  are  found 
by  tuning  for  an  inaudible  null  with  both  the  resistivity 
(i.e.  amplitude)  and  phase  dials.  The  accuracy  of  the 
instrument  depends  upon  distance  to  the  transmitter, 
as  great  distances  can  cause  much  noise  and  limit  the 
null  detection.  Since  Cutler,  Maine,  is  within  300  km 
of  the  survey  sight  and  operates  at  a  radiated  power 
of  1  MW,  noise  was  not  a  limitation. 


Airborne 

The  Barringer  Research  Ltd.  E-PHASE*  system 
was  used  for  the  aerial  survey  (Barringer  1972-73). 

*  E-PHASE  Is  the  trade  mark  of  the  Barringer  Ltd.  sys¬ 
tem. 


This  system  employs  a  Britton  Norman  Islander,  a  short 
takeoff  and  landing  (STOL)  aircraft,  specially  modified 
to  include  a  nose  stinger  on  which  horizontal  and 
vertical  dipole  antennas  are  mounted  as  shown  in 
Figure  2b.  The  E-PHASE  system  is  capable  of  moni¬ 
toring  as  many  as  four  frequencies.  For  this  survey 
only  one  frequency  was  used,  station  NAA  (17.8 
kHz),  located  near  Cutler,  Maine.  The  survey  was 
flown  with  an  average  flight  line  orientation  of  55° 

(true  azimuth)  for  an  average  coupling  angle  of  80.5° 
with  the  transmitter.  The  survey  altitude  of  the  air¬ 
craft  varied  around  1 50  m.  The  mean  flight  line 
spacing  was  0.4  km. 

Airphoto  mosaics  and/or  topographic  maps  were 
used  for  navigation.  A  flight  path  camera  obtained 
continuous  photographic  coverage  along  the  flight 
lines.  Reference  points  (fiducials)  were  placed  in  the 
flight  path  recovery  photography  and  at  correspond¬ 
ing  locations  on  data  stored  on  an  analogue  (pen  trace) 
recorder.  Manually  triggered  fiducials  were  used  to 
note  points  on  the  ground,  such  as  streams  and  lakes, 
as  well  as  the  start  and  end  of  flight  lines.  Each  day 
the  flight  path  recovery  film  was  developed,  and  re¬ 
corded  data  were  printed  and  inspected  by  Barringer 
Research  Ltd.  to  ensure  that  all  information  was  ob¬ 
tained.  In  the  case  of  missing  or  unusable  photo  cover¬ 
age,  malfunction  of  the  instrumentation,  or  transmitter 
problems,  portions  of  lines  or  entire  flights  were  flown 
again  (Barringer  1975). 

Additional  information  on  the  E-PHASE  system, 
such  as  calibration,  navigation  and  data  reduction,  has 
been  discussed  by  Hoekstra  et  al.  (1974)  and  by 
Palacky  and  Jagodits  (1975).  A  review  of  the  theory 
of  electromagnetic  resistivity  surveying  is  given  in 
Appendix  A. 

Magnetometer  survey 

A  special  tail  stinger  was  installed  on  the  survey  air¬ 
craft  to  accommodate  the  magnetometer  sensor.  The 
airborne  magnetometer,  which  records  variations  in 
the  earth’s  magnetic  field,  and  the  reference  ground 
base  unit,  measuring  the  total  field,  both  utilize  the 
principle  of  proton  precession.  A  Barringer  Research 
airborne  AM-104  proton  precession  magnetometer  was 
used  to  determine  the  variations  of  the  earth 's  magnetic 
field.  The  resolution  of  this  instrument  is  1  gamma,* 
with  a  cycling  rate  of  1.1  s.  The  magnetic  data  were 


*  7  gamma  =  7  x  70*®  tesla. 


recorded  on  an  analogue  recorder  and  on  magnetic 
tape. 

A  GM-1 23  field  proton  precession  magnetometer, 
manufactured  by  Barringer  Research  Ltd.,  was  loca¬ 
ted  at  Madawaska,  Maine  (32  km  east  of  Ft.  Kent) 
and  used  as  the  station  magnetometer.  During 
magnetic  storm  activity,  which  occurred  between 
6  and  8  October  1975,  the  surveying  was  suspended. 

In  Appendix  B  a  brief  review  of  magnetic  survey¬ 
ing  is  given. 

RESULTS 

Ground  control  study 

During  the  airborne  survey  of  the  test  sector  ground 
control  was  established  by  measuring  the  apparent 
resistivity  and  phase  by  the  surface  impedance  method. 
Readings  were  made  every  60  m  along  the  traverses 
shown  in  Figure  3,  which  are  superimposed  upon 
the  bedrock  geology  as  mapped  by  Boone  (1 962). 

The  traverses  primarily  followed  old  logging  roads, 
except  for  line  A  which  cut  along  the  steep  ridges  of 
Gardner  Mountain.  A  few  sample  readings  were  taken 
at  the  pond  shores  where  the  apparent  resistivity  in¬ 
variably  dropped  to  less  than  200  ohm-m. 

The  First  three  columns  of  Table  I  present  a  statis¬ 
tical  summary  of  the  data  collected  over  the  various 
rock  types.  Bedrock  was  rarely  more  than  several 
meters  below  the  surface  along  all  of  these  traverses. 
The  large  dispersions  in  the  resistivities,  noted  by  the 
standard  deviations  in  parentheses,  primarily  result 
from  variations  in  bedrock  properties.  The  highest 
resistivities  were  obtained  over  the  intrusive  rocks. 
Dispersion  in  phase  was  largest  for  the  slate,  with 
readings  varying  from  3°  to  54°.  Over  the  intrusive 
rocks  the  phase  quickly  stabilized,  as  indicated  by  the 
decrease  in  standard  deviation. 

For  direct  comparison  of  the  results  of  the  ground 
and  airborne  surveys,  the  ground  data  must  be  modi¬ 
fied,  based  on  the  same  assumptions  used  in  calcu¬ 
lating  the  airborne  values.  The  ground  data  are  re¬ 
computed  using  only  the  quadrature  component  of 
the  surface  impedance  and  an  assumed  phase  angle  of 
45°  (eq  A2,  App.  A).  A  normalized  distribution  of  the 
modified  ground  data  is  plotted  in  Figure  4,  showing  a 
clear  distinction  between  the  slate  and  intrusive  rocks 
in  the  control  area.  (A  summary  of  the  modified 
data  is  provided  in  the  last  column  of  Table  I).  The 


Figure  3.  Ground  control  traverses  along  which  surface  Im¬ 
pedance  measurements  were  made,  shown  superimposed  upon 
the  bedrock  geology  as  mapped  by  Boone  ( 1 962). 


Table  I.  Apparent  resistivity  data  obtained  during  ground  survey. 
The  last  column  presents  modified  ground  resistivity  values,  based 
on  the  assumptions  used  in  calculating  airborne  values.  Standard 
deviations  for  the  data  are  shown  in  parentheses. 


Bedrock  type 

Number  of 
samples 

Mean  apparent 
resistivity 
(Sl-m) 

Mean 

phase 

(°J 

Mean  quadrature 
resistivity 
(SUn) 

Seboomook  slate 

94 

5481  (0.89) 

23.93  (0.43) 

1871  (0.85) 

Granodiorite  Section  A 

22 

9009  (0.44) 

28.86  (0.15) 

4374  (0.47) 

Granodiorite  Section  B 

35 

7086  (0.70) 

35.90  (0.17) 

4597  (0.62) 

Syenite 

8 

5937  (0.56) 

43.25  (0.11) 

5342  (0.53) 
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Figure  4.  Normalized  distribution  of  the  ground  data  showing  a  dear  dis¬ 
tinction  between  the  slate  and  intrusive  rocks. 


Figure  5.  Airborne  resistivity  data  for  the  various  mapped  rock  types  (Dss- 
cyclically  bedded  gray  slate  and  sandstone,  Ds-gray  slate  and  minor  graywacke). 
Values  were  obtained  by  sampling  at  regular  intervals  along  the  flight  lines. 


suppression  of  values  for  the  slate  is  due  to  the  low 
range  of  phase  values  caused  by  the  till  cover  and/or 
the  extremely  high  permittivity  of  slate  (Parkhomenko 
1967). 

When  resistivity  values  in  Figure  4  and  the  last  column 
of  Table  I  are  compared  to  those  obtained  along  the 
flight  lines  of  the  airborne  survey  shown  in  Figure  5, 
the  agreement  may  be  seen  between  the  ground  and 
airborne  observations.  Slate  values  in  Figure  4  should 
be  compared  to  values  for  the  same  rock  type,  shown 
as  Ds  values  in  Figure  5. 

The  rapidly  changing  topography  in  the  whole  sur¬ 
vey  area  made  it  unfeasible  to  always  maintain  a  con¬ 
stant  altitude.  Therefore,  as  mentioned  previously, 
a  multi-elevation  test  was  conducted  in  a  small  test 
sector  to  determine  the  effect  of  survey  altitude, 


necessitated  by  topographic  changes,  upon  resolution. 
Two  separate  surveys  were  flown  at  mean  altitudes  of 
1 50  m  and  300  m,  both  surveys  containing  the  same 
number  of  flight  lines. 

In  Figure  6  the  comparisons  between  the  contoured 
data  from  both  surveys  can  be  made.  The  general  dis¬ 
tribution  of  resistivity  anomalies  remains  unchanged 
while  the  detail  has  decreased  at  the  higher  altitude. 

At  this  higher  altitude  there  is  a  greater  emphasis  of  the 
4000  ohm-m  contour  over  the  intrusives,  which  is  in  near 
agreement  with  the  mean  values  in  the  last  column  of 
Table  I  and  is  indicative  of  the  averaging  effect  of  altitude. 

In  general  an  excellent  level  of  correlation  was 
found  between  the  airborne  and  ground  readings  for 
the  test  sector,  ensuring  the  quality  of  the  remaining 
resistivity  survey. 
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No  ground  studies  were  conducted  for  the  airborne 
magnetometer  survey,  since  past  research  and  employ¬ 
ment  of  this  technique  have  already  established  its 
value. 

VLF  survey 

The  most  noticeable  features  of  the  VLF  data  are 
the  numerous  tightly  clustered,  closed-contour  re¬ 
sistivity  highs  shown  in  Figure  7.  This  pattern  reflects 
the  complicated  surficial  and  bedrock  geology  of  the 
area.  In  general  the  highest  resistivity  areas  are  associ¬ 
ated  with  areas  of  major  relief.  This  correlation  is 
based  on  the  fact  that  these  are  the  locations  where 
the  conductive  til!  is  the  thinnest  and  values  reflect 
the  generally  more  resistive  bedrock  which  is  near 
the  surface.  A  generalized  topographic  map  is  pro¬ 
vided  (Fig.  8)  as  an  aid  in  illustrating  that  the  bedrock 
types  can  best  be  examined  in  areas  of  greatest  relief. 
Virtually  all  the  areas  of  high  resistivity  (>  3000  ohm- 
m)  agree  with  units  mapped  as  bedrock,  or  bedrock 
with  a  till  cover,  in  the  earlier  surficial  geologic  study 
conducted  by  CRREL  (Mckim  and  Merry  1975). 

A  quantitative  demonstration  of  these  observations 
is  shown  in  Table  II.  The  table  gives  the  distribution 
of  areas  with  greater  than  3000  ohm-m  resistivity  in 
relation  to  the  topographic  features  they  contain. 
Seventy-three  areas  were  delineated,  none  containing 
only  a  valley.  Approximately  70%  of  the  areas  con¬ 
tained  only  flanks  and/or  ridges. 

Table  II.  VLF  apparent  resistivity  contours  greater 
than  3000  ohm-m. 


Topographic  situation 

Number 

Percentage 

Ridge 

14 

19.2 

Flank 

17 

23.3 

Ridge  and  flank 

20 

27.4 

Ridge,  flank  and  valley 

16 

21.9 

Flank  and  valley 

6 

8.2 

Valley 

0 

0 

Totals  73  100.0 

It  is  also  apparent  from  a  comparison  of  Figure  7 
and  Figure  8  that  many  mountainous  areas  are  also  of 
relatively  low  resistivity  (<  3000  ohm-m).  Therefore, 
the  lack  of  conductive  overburden  is  not  a  major  in¬ 
fluence  upon  resistivity  distinctions  within  these  areas. 


The  most  plausible  explanations  for  the  resistivity 
distinctions  among  mountain  types  are  the  following: 

1)  There  is  a  difference  in  bedrock  type. 

2)  Certain  flights  passed  directly  over  mountain 
peaks  where  it  is  known  (Harrison  et  al.  1971 )  that 
the  vertical  electric  field  increases  because  of  the 
sharp  features,  thereby  suppressing  resistivity  values. 

The  second  explanation  is  applicable  to  individual 
peaks  such  as  in  the  Hafey,  Gardner,  and  Deboullie 
mountains.  Along  these  high  elevations,  ground 
readings  (granodiorite  Section  A  in  Table  I  and 
Sellmann  et  al.  1975)  that  were  modified  for  compari¬ 
son  to  the  airborne  method  indicated  a  suppression 
of  airborne  values.  However,  most  flightlines  in 
mountainous  areas  passed  over  the  flanks,  rather  than 
peaks. 

Bedrock  geology  and  resistivity 

When  the  reconnaissance  bedrock  geology  of  this 
area  (Boudette  et  al.  1966  and  Boone  1962)  is  ex¬ 
amined,  general  distribution  patterns  become  apparent. 
As  shown  in  Figure  9,  in  the  central  part  of  the  study 
area  the  bedrock  units  occur  in  broad  bands,  while  to 
the  northern  end  of  the  area  these  bands  become  more 
contorted  and  narrow,  including  a  greater  contrast  in 
rock  types.  A  large  intrusive  rock  unit  occurs  near  the 
southern  limit  of  the  area. 

The  southern  and  central  parts  of  the  region  were 
characterized  by  three  main  units  of  lower  Devonian 
age.  The  southernmost  unit,  Ds,  was  defined  as  gray 
slate  and  minor  graywacke.  The  southern  limit  of  this 
slate  unit  is  not  mapped,  but  it  probably  extends  be¬ 
yond  the  limits  of  the  bedrock  reconnaissance  study, 
to  the  southern  extent  of  the  present  study.  The 
central  part  of  the  area  was  mapped  as  a  cyclically 
bedded  gray  slate  and  sandstone  (Dss).  North  of  this 
unit  another  lower  Devonian  unit  (Dsg),  composed 
of  graywacke  and  gray  slate,  was  defined.  North  of 
these  large  units  the  structure  becomes  more  com¬ 
plex  and  rock  units  range  in  age  from  Ordovician  to 
Devonian,  A  greater  range  in  rock  types  also  occurs, 
including  limestone,  orthoquartzite,  and  siltstone. 

The  intrusive  body  found  in  the  southern  part  of  the 
area  is  composed  of  granodiorite  and  syenite. 

When  the  resistivity  contours  are  placed  upon  Figure 
9  several  correlations  become  apparent.  The  highest 
resistivity  values  (>  6000)  are  associated  with  the  intru¬ 
sive  rocks  in  the  Gardner  Mountain  area.  The  next 
most  obvious  correlation  is  seen  when  areas  of  greater 
than  3000  ohm-m  are  examined.  A  belt  containing 


Figure  8.  Generalized  topographic  map  with  selected  contour  Intervals  (ft)  and  drainage 
patterns  (after  USGS  topographic  maps,  1:62500  series). 
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Figure  9.  Bedrock  distribution  map  (after  Boudette  et  at.  1966,  Boone  1962). 
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many  high  resistivity  zones  is  noticeable  in  the  central 
part  of  the  study  area,  coinciding  with  the  mapped 
limits  of  the  Dss  bedrock  unit.  In  this  belt,  the  high 
resistivities  are  found  in  areas  of  high  relief.  The 
values  commonly  exceed  3000  ohm-m  and  locally 
reach  a  maximum  of  5500  ohm-m. 

The  remaining  rock  types,  such  as  the  Ds  geologic 
unit,  have  very  low  resistivities,  with  few  areas  having 
values  that  exceed  3000  ohm-m.  This  occurs  despite 
the  wide  range  of  rock  types  and  large  number  of 
areas  of  high  relief  where  bedrock  is  near  the  surface. 
The  resistivity  data  from  the  Dss,  Ds,  and  intrusive 
units  in  Figure  5  show  these  same  trends.  The  intru¬ 
sive  unit  with  the  highest  resistivities,  Dss,  has  more 
than  20%  of  its  values  greater  than  3000  ohm-m  while 
less  than  2%  of  the  values  in  the  Ds  unit  exceed  3000 
ohm-m. 

Aeromagnetic  survey  results 

The  aeromagnetic  survey  did  not  outline  anomalies 
usually  associated  with  abnormally  high  magnetic 
mineralization  in  either  the  main  study  area  or  along 
flight  line  pairs  flown  along  the  major  drainage  net¬ 
works.  In  general,  the  magnetic  intensity  data  were 
featureless,  with  the  only  significant  anomaly  being 
associated  with  the  known  intrusive  rocks  in  the 
southern  part  of  the  area.  This  anomaly  is  shown  in 
Figure  10  for  comparison  to  the  resistivity  data  and 
the  bedrock  geology.  Fifty-seven  thousand  gammas 
have  been  subtracted  from  all  values  to  eliminate  the 
background  intensity  of  the  earth.  Modal  analysis  of 
the  granodiorite  and  syenite  (Boone  1962)  indicated 
that  magnetite  is  a  common  and,  in  some  cases, 
abundant  accessory  mineral.  This  fact  alone  can 
account  for  this  distinctive  feature,  although  the 
patterns  and  intensities  seen  are  considered  normal 
for  these  rock  types. 

CONCLUSIONS 

The  combined  results  of  the  aeromagnetic  and 
E-PHASE  survey  allow  several  conclusions  to  be  drawn 
concerning  sources  of  rock  suitable  for  construction 
purposes.  These  conclusions  are  based  on  the  fact 
that  the  most  suitable  rocks  will  be  the  most  electri¬ 
cally  resistive,  with  further  distinctions  based  on 
mineralization  to  be  inferred  from  the  presence  of 
magnetic  data. 


The  conclusions  are: 

1 .  The  highest  resistivity  values  are  associated 
with  the  known  Gardner  Mountain  intrusives.  This 
is  also  the  location  of  the  only  significant  magnetic 
anomaly.  The  correlation  of  the  VLF  and  magnetic 
data  in  this  area  was  shown  in  Figure  10. 

2.  The  lack  of  similar  resistivity-magnetic  correla¬ 
tions  in  other  locations  suggests  little  or  no  other 
occurrence  of  these  rock  types. 

3.  Local  resistivity  highs  in  the  Dss  unit  suggest 
locations  for  field  study  as  possible  sites  of  rock  types 
suitable  for  construction  purposes.  The  Boudette  et 
al.  (1966)  study  indicates  this  as  an  area  of  slate  and 
sandstone.  The  distribution  of  these  rock  types  is 
not  known;  therefore,  it  can  only  be  assumed  that 
the  sandstones  may  be  the  most  resistive  rock  type. 

No  magnetic  anomalies  occur  in  this  area. 
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APPENDIX  A.  THEORY  OF  ELECTROMAGNETIC 
RESISTIVITY  SURVEYING 


RESISTIVITY  SURVEYING 

The  physical  property,  electrical  resistivity,  spans 
a  large  range  of  magnitude.  Values  ranging  from  10*7 
ohm-m  for  metallic  elements  to  1010  ohm-m  for  some 
granitic  rocks  are  well  documented  (Parkhomenko 
1967).  The  most  common  range  encountered  in  d.c. 
and  electromagnetic  surveying  is  approximately  10  to 
10,000  ohm-m.  Within  this  range  fall  many  materials 
of  economic  and  geologic  importance  that  can  be  de¬ 
tected  as  a  result  qf  contrasts  in  resistivity  with  ad¬ 
jacent  materials. 

A  primarv  influence  on  the  resistivity  of  sediments 
is  the  amount  of  ions  associated  with  adsorbed  surface 
water  on  particles.  Therefore,  there  is  a  correlation 
between  resistivity  and  grain  size.  In  some  earth 
materials,  such  as  crystalline  rocks,  a  more  important 
factor  may  be  the  amount  of  ions  dissolved  within  the 
pore  water.  The  general  resistivity  range  for  most 
earth  materials  is  illustrated  in  Figure  A1  (Culley  et  al. 
1975,  Hoekstraand  Delaney  1973).  Other  factors, 
such  as  temperature,  pressure  and  ice  content,  also 
influence  resistivity  but  are  not  significant  factors  in 
the  present  study. 

Radiowaves  propagating  over  the  earth’s  surface  are 
influenced  by  the  electrical  properties  of  the  subsurface 
materials.  By  comparing  the  various  field  components 
of  the  radio  surface  wave  it  becomes  possible  to  obtain 
ground  resistivity  data. 

Since  literature  on  the  interaction  between  radio¬ 
waves  and  the  earth  is  available  (Wait  1962,  Frisch- 
knecht  1973,  Eliassen  1956),  only  a  short  discussion 
of  the  principles  is  presented  here. 

The  electromagnetic  ground  wave  field  vectors  in 
the  far-field  of  a  vertically  polarized  transmitter  are 
shown  in  Figure  A2.  At  the  ground  surface  there  are 
three  field  vectors:  the  horizontal,  radially  oriented 


electric  field  £x,  the  horizontal,  azimuthzlly  oriented 
magnetic  field  Hv,  and  the  vertical  electrical  field  £z 
{x,y,  z  refer  to  a  local  Cartesian  coordinate  system). 

All  three  field  vectors  decay  equally  in  amplitude  with 
increasing  distance  from  the  transmitter.  In  the  ground 
wave  mode  the  relative  amplitudes  and  phases  of  these 
fields  at  VLF  are  not  influenced  by  the  total  path  of 
propagation  but  only  by  local  resistivity  conditions. 

At  distances  approaching  1000  km  sky  mode  trans¬ 
mission  becomes  important.  (For  the  survey  discussed 
in  the  text,  the  proximity  of  the  test  site  to  the  t>  ans- 
mitter  (300  km)  ensured  that  only  the  ground  wave 
was  present  and.hence  that  there  was  no  possibility 
of  ionospheric  interference.) 

Because  of  the  large  refractive  index  of  the  ground 
at  radio  frequencies,  a  near  vertical  wave  propagates 
into  the  ground  with  horizontal  vectors  Wy  and  £x. 

In  the  ground  the  amplitudes  of  £x  and  Wy  attenuate 
with  depth,  and  in  homogeneous  ground  the  distance 
over  which  the  field  decays  to  37%  of  its  surface 
value  is  called  the  skin  depth  of  the  radiation.  Figure 
A3  illustrates  the  dependency  of  the  skin  depth  on 
the  frequency  of  the  transmitted  signal  and  the 
resistivity  encountered  for  a  homogeneous  earth. 

The  skin  depth  of  the  radiation  is  an  important 
parameter,  since  it  approximately  indicates  the  depth 
to  which  information  is  obtained  by  a  measurement 
at  the  surface. 

Both  £x  and  Hy  are  continuous  across  the  ground 
surface,  but  £z  becomes  negligible  in  the  ground. 

£x  can  therefore  be  measured  by  a  dipole  antenna 
above  the  earth’s  surface  or  in  the  ground  by  measur¬ 
ing  the  field  strength  between  two  probes.  //y  is 
measured  by  a  coil  located  near  the  surface,  and  £z 
is  measured  with  a  vertical  dipole  antenna  above  the 
ground. 

The  basis  for  obtaining  a  local  measurement  of 
ground  resistivity  is  illustrated  in  Figure  A4,  where 
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After  Cut  ley  et  at.  (1975). 


From  Hoekstra  and  Delaney  ( 1973). 


Figure  A 1.  Examples  of  resistivity  ranges  for  common  earth  materials.  Soil  types  are  described  using  the  Unified 
Soil  Classification  System. 


Patriot# 

Transmitter 


Figure  A2.  Electromagnetic  field  components  of  a  vertically  polarized  radio  sur¬ 
face  wave. 
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Figure  A3.  Skin  depth  of  electromagnetic  plane  waves  as  a  function  of  fre¬ 
quency  with  varying  ground  resistivity  p. 


E -  PHASE  RESPONSE 


Figure  A4.  Schematic  representation  of  the  local  changes  In  the 
magnitude  of  the  electromagnetic  field  vectors  of  a  vertically 
polarized  groundwave  when  earth  materials  having  different  re¬ 
sistivities  are  traversed. 


a  wave  propagates  over  a  change  in  ground  conditions. 
Changes  in  local  subsurface  conditions  cause  perturba¬ 
tions  in  amplitude  and  phase  of  £„,  while  local  changes 
marginally  affect  Ez  and  Hy.  Therefore,  by  measuring 
the  ratio  of  EJHy  (called  surface  impedance)  or  EJEt 
(called  wave  tilt)  a  measurement  is  obtained  of  the 
local  electrical  resistivity. 

The  computation  of  ground  resistivity,  at  present, 
neglects  the  effects  of  the  permittivity  of  the  ground, 
and  there  are  circumstances  when  such  procedures  may 
introduce  substantial  errors  (Olhoeft  1975).  However, 


in  many  applications  the  delineation  of  changes  in 
subsurface  conditions  is  the  primary  object,  and  a 
determination  of  the  exact  value  of  ground  resistivity 
is  less  important. 

Over  the  past  decade  several  schemes  for  measuring 
the  resistivity  of  ground  from  radiowaves  by  using  the 
principles  discussed  above  have  been  tried.  Two 
approaches  have  proven  to  be  practical:  a  ground- 
based  technique  (Collett  and  Becker  1967),  and  an 
airborne  technique  (Barringer  1972, 1973).  The 
ground  technique  relies  on  a  measurement  of  surface 
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impedance  (£x///y),  while  the  airborne  method  uses 
the  wave  tilt  phenomenon  (£x/£z)  for  plane  waves  at 
grazing  incidence.  For  plane  waves  at  grazing  incidence 
over  a  smoothly  layered  earth,  the  two  measurements 
are  related  by  the  free  space  impedance,  377  ohms,  as 
follows: 

W  =  Z,/377 

where  W  is  the  wave  tilt  and  Zf  the  surface  impedance. 

Since  Z,  and  W  are  complex  quantities,  they  both 
have  an  associated  amplitude  and  phase.  The  ampli¬ 
tude  is  converted  to  "apparent”  resistivity  by  using 
the  formula  for  a  homogeneous  earth  (Wait  1962) 
such  that 

Pgpp  =  Ef/2nf(i0  =  W7l2irfe0  (A1) 

where  papp  is  the  apparent  resistivity  measured  in  ohm- 
meters,  f  the  operating  frequency  in  Hz,  p0  the  mag¬ 
netic  permeability  of  free  space  and  eQ  the  free  space 
permittivity.  Theoretical  modifications  in  apparent 
resistivity  and  phase  due  to  layering  of  materials  with 
differing  resistivities  are  well  documented  (Wait  1962). 
In  general  when  resistivity  increases  with  depth,  phases 
are  below  45°,  and  when  resistivity  decreases  with  depth 
the  phases  are  greater  than  45°.  Phases  of  45°  are  indi¬ 
cative  of  uniform  resistivity  with  depth. 

For  the  airborne  system  only  the  quadrature  value 
of  wave  tilt  can  be  measured,  because  aircraft  roll 
instability  causes  an  in-phase  coupling  of  the  vertical 
field  with  the  (offset)  horizontal  antenna.  Therefore, 
present  data  processing  assumes  an  arbitrary  phase 
angle  of  45°  between  Ex  and  £z  and  the  computation 
of  apparent  resistivity  uses  the  formula: 

P,  Pp  =  W2lnfe0  (A2) 

where  W  =  quad  (£x/£z). 

This  inability  to  distinguish  the  separate  effects  of 
phase  and  amplitude  limits  the  degree  of  airborne  sub¬ 
surface  interpretation.  Therefore,  it  often  becomes 
necessary  to  conduct  preliminary  ground  studies  to 
determine  if  the  combined  effect  of  phase  and  ampli¬ 
tude  will  still  allow  extensive  airborne  differentiation 
of  material  type. 

Ground  measurements  are  performed  with  the  sur¬ 
face  impedance  (EJHy)  technique.  The  technique 
must  be  used  for  groundbased  surveys  because  Ez  is 


disturbed  by  vegetation  and  is  therefore  an  unre¬ 
liable  reference.  In  the  airborne  survey  it  becomes 
necessary  to  use  the  quadrature  wave  tilt  method 
to  eliminate  the  undesirable  coupling  mentioned 
above.  In  this  case  £z  is  more  reliable  since  the  ob¬ 
servation  is  made  well  above  the  ground  surface,  sub¬ 
stantially  removed  from  vegetation  effects.  However, 
mountain  ridges  are  known  to  sometimes  enhance 
this  reference  field  (Harrison  etal.  1971)  thereby 
depressing  the  actual  resistivity  values  but  not  alter¬ 
ing  regional  resistivity  patterns,  as  is  discussed  in 
Results. 

Since  the  radiation  wavelength  used  is  so  great 
(>  1 5  km  at  VLF),  a  flight  altitude  of  1 50  m  is 
equivalent  to  less  than  0.01  wavelength.  However, 
the  resolution  in  airborne  mapping  of  resistivity 
anomalies  is  reduced,  compared  to  ground  measure¬ 
ments. 


APPENDIX  B.  MAGNETIC  SURVEYING 


Magnetic  data  are  a  result  of  both  the  intensity  of 
the  earth’s  magnetic  field  and  the  magnetic  properties 
of  minerals  present  in  the  surveyed  area.  Magnetic 
anomalies  are  disturbances  in  the  earth’s  magnetic 
field  that  are  caused  by  changes  in  the  amount  and 
type  of  magnetic  minerals  associated  with  various  rock 
types.  The  magnetic  minerals  usually  responsible  for 
these  disturbances  are  magnetite,  pyrrhotite,  ilmenite, 
and  hematite. 

Magnetic  disturbances  present  themselves  in  various 
ways,  such  as  changes  in  field  direction  (usually  termed 
a  "dip"  or  change  in  declination),  and  as  changes  in 
total  field  strength  or  intensity.  The  most  common 
method  of  presentation  is  contouring  lines  of  total 
magnetic  intensity  above  a  background  value  measured 
at  some  reference  station.  Values  are  presented  in 
terms  of  gammas  above  the  reference  station  value, 
one  gamma  being  equal  to  10'5  gauss.  The  earth’s 
magnetic  field  varies  with  position  and  time  of  year 
but  generally  falls  in  the  magnitude  range  of  0.5  to  1 
gauss.  Contours  do  not  necessarily  outline  a  rock 
formation  but  can  be  highly  indicative  of  the  rock 
types  present. 

Additional  information  on  magnetometer  survey¬ 
ing  can  be  found  in  the  texts  by  Grant  and  West  (1 965) 
and  Ward  (1967). 


